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ABSTRACT 



Context. Observations indicate a strong link between star formation and black hole (BH) growth, but some questions remain unan- 
swered: whether both activities are coeval or whether one precedes the other, what their characteristic timescales are, and what kinds 
of physical processes are responsible for this interplay. 

Aims. We examine stellar and BH masses (M, and M s #) in z~2 active systems at the peak of their AGN or star formation activity 
to investigate how they are linked and whether AGN radiative or else radio power provides a feedback mechanism that regulates the 
stellar growth in these systems. 

Methods. We analyze the infrared (IR) spectral energy distributions of radio, sub-millimeter and mid-IR selected AGNs at z ~1— 3 
and constrain their stellar and AGN luminosities using AGN and host-galaxy templates. 

Results. We find evidence of increasing stellar light, thereby decreasing the AGN mid-IR power going from mid-IR selected AGNs, 
to radio galaxies, and to sub-millimeter AGNs. This trend can be explained by either decreasing Eddington ratios or increasing offsets 
from the local M BH -M, relation. All systems are characterized by high star formation rates regardless of their different AGN powers, 
thus neither AGN radiative power nor AGN-driven radio activity seems to influence the star formation rate in the selected AGNs. We 
discuss two possible evolutionary scenarios that might link these three AGN classes. 

Key words. Galaxies: active - Galaxies: evolution - Galaxies: high-redshift - quasars - Infrared: galaxies 



1. Introduction 

A new paradigm has recently emerged in which galaxies and 
their black holes (BHs) grow in tandem, and BHs are respon- 
sible for regulating star formation (SF) through feedback pro- 
cesses and thus shaping many of the basic rel a tionships between 
Crotonetajl l2Q06h iMerloni & Heinzl 
Lin et all l2003t ISomerviUd I2004I). 



galaxy properties (e.g 



y p 

I2007t iBalogh et alJ 12001 



Feedback processes are thought to be driven by AGNs, and they 
manifest themselves in outflowing winds or radio jets that can 
heat a significant fraction of the gas even, from the most mas- 
■ sive dark-matter halos, and terminate the S F process in galax- 
ies dDi Matteo et alj2005tlCroton et al.ll2006l) . 

Current evolutionary models predict that AGN- 
feedbac k takes place in o bscured and extremely lu minous 
AGN s dHopkins et all 120051) or in radio-loud AGNs dCrotonl 
2006), especially at z >2 where both BH accretion and SF 
activity peak (e.g. iGiavalisco et all 12004 iMarconi et all 12004). 
Although not in large numbers, these kinds of objects have 
been found in wide-area radio and infrared surveys, where the 
effects of obscuration are minimized and AGN signatures are 
clearly identified. Moreover, thanks to the advent of the Spitzer 
Space Telescope, it has become possible to study the multi- 



wave l ength properties of these objects in deta il (Ogle et al 



2008). 



2006; ISevmour et ail 120071; ISajina et alJ |2007b; Pollett a et al 



In this Letter, we compare the relative luminosity produced 
by stars and by AGN-heated dust in three classes of AGNs, high- 
z radio galaxies (HzRGs hereinafter), sub-millimeter-detected 
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AGNs (SMG/AGNs hereinafter), and IR-selected obscured 
QSOs (IsOQs hereinafter), to investigate whether these three 
classes represent members of the same population caught dur- 
ing different evolutionary stages, or whether they are the result 
of different histories and physical conditions. 



2. Obscured AGNs at z ~2 

HzRGs have been discovered due to to their large radio 
fluxes (>10mJy) in shallow wide surveys ( e.g. NVSS, WENSS 
ICondon et all 1 19981: iRengelink et alj|1997l) . Since they are ex- 
tremely rare, it is necessary to sample large volumes to find 
them. They are usually characterized by high excitation emis- 
sion lines in the optical and strong mid-IR (MIR) emission, and 
their radio emission is dominated by non-thermal radiation from 
radio lobes. In the X-ray, they show both the contribution from a 
jet and an ofte n absorbed component associated with the accre- 
tion disk (e.g. lSevmour et alJl2Tj07l:lHardcastle et al.ll2006l) . 

IsOQs have been been discovered as bright (>lmJy) 
24pm sources in Spitze r surveys with fa i nt optical coun - 
terparts (F Um /F R >200; iHouck et all 12004 lYan et all 120051) . 
Their MIR SEDs are usually consistent with power-law 
or convex spectral shapes rapidly rising towards longer 
wavelengt hs, and they often displ a y absorption due to 
silicates ( Alonso-Herrero et al J [2006; Weedma n et alJ 2006; 
iPolletta et al.ll2006Ll2008h ~ 

SMG/AGNs were first discovered in sub-millimeter sur- 
veys, and their AGN nature was revealed by their X-ray emis- 
sion (Alexander et al.1 1200 5b). X-ray observations indicate that 
these galaxies host moderately luminous and heavily absorbed 
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AGNs. Their optical-IR emission is mainly dominated by the 
host-g alaxy, and they are characterized by inten se starburst ac- 
tivity dBorvs et al]|2005b IXTexander et alll2005blla|). 

Fo r this work we selected 17 HzRGs from Seymour et al.l 
d2007l) . 1 1 X-ray d etected SMG/AGNs from lBorvs et al l(l2005l) . 
and 21 IsOQs from lPolletta et al.l d2008l) . The three samples only 
include sources at l<z<3. The median and redshift range of 
each sample are reported in Table Q] All luminosities were de- 
rived assuming Hq = 70km s _1 Mpc~', Qm=0.3, and Qa=0.7. 
Throughout the paper, the uncertainties associated to the median 
values correspond to the average deviation from the median. 

2. 1 . Host stellar light and AGN accretion power 

A common problem in AGN studies is to estimate the host- 
galaxy contribution to the bolometric luminosity and, especially, 
the stellar mass. Since the stellar emission peaks in the near-IR 
(NIR), typically at 1.6 pm in the rest frame (or H-band), and the 
ratio between NIR luminosity and stellar mass (M.) are char- 
acterized by little dispersion, the NIR luminosity is a proxy of 
the bulge mass or the luminosity of the host-galaxy. Note that 
the L(NIR)/M» depends on the SF history and age of the sys- 
tem, e.g. younger stellar populations can yield 30 % lower stellar 
mass estimates than older stellar populations (e.g. lSevmour et all 
|2007|) . Since all selected objects are at similar redshifts and their 
AGN component does not dominate at NIR wavelengths, rea- 
sonable estimates of stellar emission and mass can be obtained 
from L(NIR). In IsOQs, the NIR stellar light is estimated af- 
ter s ubtracting the AGN c omponent. We used the NIR estimates 
from lPolletta et all d2008l) . where the IR SEDs are modeled with 
a combination of ho st and torus models. Fo r the HzRGs, we used 
the estimates from Sey mour et al.l d2007l) where the contribu- 
tion from AGN-heated dust to the NIR is also removed. NIR 
lumin osities were deriv ed for the SMG/AGNs by fitting their 
SED (Bor vs et al.l 12005) with galaxy templates. Stellar masses 
were d erived from L(H) following the method in lSevmour et alj 
(2007). The estimated NIR luminosities and corresponding stel- 
lar masses are reported in Figure Q] 

The IR SEDs of HzRGs and IsOQs indicate that their MIR 
emission is dominated by AGN-heated dust. Since the MIR (i.e., 
at 5 pm in the rest frame) is not significantly affected by dust ex- 
tinction, and this is where the AGN-heated dust radiation peaks 
and where the emission from star forming galaxies shows a min- 
imum, it can be considered a good proxy for AGN radiative 
power. For SMGs, the MIR luminosity is contaminated by the 
host-galaxy light, therefore it only provides an upper limit to 
the AGN power. M IR luminosities are available for the HzRGs 
from the literature (Sevmour et al. 2007). For the IsOQs, they 
are estimated using the models in Polle tta et al. (2008). For the 
SMG/AGNs, they are derived by fitti ng their SEDs inclu ding 
24jum data, with galaxy templates (see Polletta et al. 2 0071) . 

The stellar NIR luminosities of all selected sources are com- 
pared with the MIR luminosities in Figure Q]Note that the lack 
of IsOQs at L (MIR)< 1O 12 L is due to the way these objects 
were selected dPolletta et al.ll2008l) . 

The selected AGNs are characterized by a broad range of 
NIR/MIR (stellar/AGN) luminosity ratios. Since the NIR lu- 
minosity is a proxy for stellar mass and the MIR luminosity 
a proxy for AGN power, the broad range of L(NIR)/L(MIR) 
(stellar/AGN) luminosity ratios indicate that the three classes 
of selected AGNs emit at different Eddington ratios or fol- 
low different M^-M, relations. These two possible interpre- 
tations are explained below and illustrated in Figure Q] We 
consider the local correlation between the bulge mass and BH 



Log(M B H/M )-Log(A)-Loq(<5) 
7.0 7.5 8.0 8.5 9.0 9.5 



1 2.5 



1 2.0 



1 1.5 



q 1 1 .0 



10.5 



10.0 



• HzRG (L H '°') 
"O HzRG (L H "° r ) 
XlsOQ (L„ s,or ) 

ASMG/AGN 




12.5 



12.0 



11.5- 



1 1 .0 ■ 



10.5 



10.0 



10.5 11.0 11.5 12.0 12.5 13.0 
Log(L(5 M m)/L ) 

Fig. 1. Comparison between the H-band (1.6 pm) and 
the MI R (5 Ami) luminositie s of three classes of AGNs, 
HzRGs (Sev mour et al. 2007]) {fu ll and open blue circles), 



SMG/A GNs dBorvs et all l2005l) {green triangles), and 



IsOQs dPolletta et alj 120081) {red crosses). The H-band lu- 
minosity of the HzRGs refers to the total {full blue circles) 
or stellar-only (open blue circles) luminosity. Downward 
pointing arrows represent upper limits to the H-band lumi- 
nosity, and left pointing arrows indicate upper limits to the 
MIR luminosity. The diagonal lines represent the expected 
luminosities for BHs of various masses, as indicated on the 
top axis, assu ming Lfff^/T,F^=/L where is derived 

from L( 5 urn) dPolletta et all 120081) . and the local M BH -M* 
relation ( Marconi & Hunt! 120031) where M, is derived from 
L(H) dSevmour et al.l 120071) (see right vertical axis) with an 
offset 6 so that higher or lower BH masses are derived for, 
respectively, S> 1 or < 1 . The names of radio-detected IsOQs are 
annotated just to the right of their location (see TableQ]). 



mass ( Marcon i & Huntl 120031) and a possible offset given by 
6, i.e. Log(M B// ) = 8.28+0.96x(Log(M«)-10.9)+Log(<5), where 
Mbh is the BH mass in M , and M» the bulge mass in M Q 
derived from L(H). The BH mass can be derived from the 
MIR luminosity assuming an Eddington ratio (L^ w fl^Edd=X) 
and a relationship between L(MIR) and the AGN bolometric 
luminosity, Log(L^' v ). The latter can be estimated using the 
following relation, Log(L^ A ') = Log(L(5/zm))+0.36±0.07, de- 
rived from the reprocess ed thermal emission produced by the 
torus (Pol letta et alj2008l) . This relation might underestimate the 
AGN bolometric luminosity because it does not include the AGN 
absorbed radiation. For compari s on, wh en using the unobscured 
AGN template from lElvis et all dl994|) , which includes optical 
and far-IR components, Log(L^ A ') = Log(L(5yum))+1.16. 

Assuming the BH and stellar masses estimated as described 
above, we derived the expected NIR and MIR luminosities for 
different values of Eddington ratio (A) and of the offsets from 
the local M^-M, relation (6) and report them in Figure[T] Note 
that all expected MIR luminosities wou ld be higher by a fac- 
tor 6.3 if the L^ A, -L(5 J um) relation from lElvis et all d 19941) was 
adopted. Note that these relations, as well the estimated masses, 
are characterized by large uncertainties. However, the uncertain- 
ties are smaller than the offsets observed among the three AGN 
samples, and all systematic uncertainties would shift all the val- 
ues by equal amounts leaving the observed dispersion and offsets 
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still present among the three AGN classes. Since our interpreta- 
tion is based on such offsets rather than on absolute values, our 
results are not affected by these uncertainties. 

The IsOQs show systematically higher Lf G , N fLEdd ratios, 
consistent with being close to Eddington-limited, while HzRGs 
and SMG/AGNs are characterized by progressively lower 
Eddington ratios. Instead of different Eddington ratios, the dif- 
ferent L(MIR)/L(H) luminosity ratios could be explained by dif- 
ferent Mflfl-M* relations with different offsets from the local 
one. Assuming that all the selected AGNs are Eddington-limited 
would imply a relatively small offset (0.1<<5<10) from the lo- 
cal Mflff-M* relation in IsOQs and larger offsets in HzRGs and 
SMG/AGNs (0.01<<5<0.5) and thus smaller BH masses than 
expecte d from the l ocal r elationship at parity of stellar mass. 
Indeed. iBorvs et all d2005l) find that SMG/AGNs are character- 
ized by lower BH masses than expected from the local rela- 
tionship assuming that they are Eddington-limited. Alternatively, 
they would lie on the M#//-M« relation if a lower Eddington ra- 
tio was assumed. In summary, the offset in L(H)-L(MIR) space 
between HzRGs and AGN/SMGs and the IsOQs can be at- 
tributed to a difference in L^J^/L^ ratios or to different offsets 
from the local M^n-M, relation. 

In the next section, we derive the median SED of the three 
samples and model them with starburst and obscured AGN tem- 
plates. 



2.2. Obscured AGN infrared SEDs 

In Figure we plot the median SEDs of HzRGs, SMG/AGNs, 
and IsOQs after normalizing them at 5 /mi in the rest frame 
at the median MIR luminosity of each group. The median 
Log(L(5 //m)/L ) is 12.47 ±0.23, 11.55+0.43, and 10.96+0.36 
for IsOQs, HzRGs, and SMG/AGNs, respectively. The me- 
dian and dispersion values are derived by regrouping the data 
points of the normalized SEDs in bins with a width of 0.12- 
0.2 in Log(/J) or larger for less than 4 data points, or smaller 
for more than 18 data point s. We also show th e median tem- 
plate of an obscured QSO dPolletta et al.ll2007l) . an d the SEP 
of th e prototypical starburst (SB) galaxy M82 dSilva et al.l 
1998). The SB template is normalized at the median stellar H- 
band luminosity of each group. The median Log(L(H)/L ) is 
10.95+0.32, 11.32+0.24, and 11.37+0.21 for IsOQs, HzRGs, 
and SMG/AGNs, respectively. The QSO template is added af- 
ter applying a certain degree of foreground extinction, Ay<8, 
and normalized so that the sum of the QSO and the SB tem- 
plates are consistent with the median SED of each group. We 
use the Galactic center extinction law by IChiar et ail ([2006) 
extinction curve because it extends to the MIR and has been 
shown to reproduce the IR spectra of AGNs (Sai ina et al.ll2007at 
iPolletta et al.ll2008h . Note that the adopted range of extinction 
values, Ay<8, corresponds to extreme dust opacities (up to hun- 
dreds of magnitudes) if a proper treatment of absorption and re- 
emission is considered as in radiative transfer models. Thus, the 
selected Ay range represents both the extinction that could be 
produced by a dusty galaxy and an optically thick torus. The 
sum of the extincted QSO template and of the SB template is 
also shown in Figure [2] 

Figure [2] shows progressively higher stellar luminosities 
and weaker AGN radiative power in IsOQs, HzRGs, and 
SMG/AGNs. Motivated by recent models that predict feedback 
on SF and BH gr owth from AGN-driven radio activity (e.g. 
ICroton eFa l. 2006), in the next section, we analyze whether the 
observed differences in stellar masses and AGN radiative power 
are related to the radio properties of the three AGN classes. 



2.3. Radio properties 



Radio fluxes are available for all H z RGs dSevm our et al. 
120071) and SMG/AGNs dBorvsetal.1 120051) and for 6 



IsOQs dBecker et alj Il995t IPolletta et al.l 120061 ISajina et all 
2007b). The radio luminosities and MIR/radio luminosity ra- 
tios of the 6 IsOQs are listed in Table 1, and their location 
shown in Figure [T] We also list median values and ranges for 
all 3 classes. For the HzRGs, we converted the radio lumi- 
nosities at 3 GHz to 1.4 G Hz assuming a pow er-law model 
with spectral index a — -0.8 dSevmouf et al. 2007), where S v oc 
V. It is well-established that the radio emission in HzRGs 
is AGN-powe red, while it is ma inly due to SF activity in 
SMG/AGNs dChapman et al.ll2004l) . The two classes show very 
different radio fluxes, luminosities and L(MIR)/L(Radio) lumi- 
nosity ratios. The HzRG L(MIR)/L(Radio) luminosity ratio is 
>10, and it is ~10 4 in SMG/AGNs (see Table 1). Based on the 
ratio between the radio and the 24/im luminosity dPonlev et al.l 
2005), we find that one source (LH_A8) is radio-loud and 
2 are radio intermediate (N2_09, and MIPS22204). For 3 of 
the 6 IsOQs with known radio fluxes, SFRs are also avail- 
able from the measured far-IR lu minosity, LH_A8, N2_09, and 
MIPS22204 dPolletta etal. 1 l2008T). In al l 3 ca ses, the SFR de- 
rived from the radio luminosity dCondonll 19921; rKennicud fl998) 
is a factor of 10 or 100 higher than derived from the FIR lumi- 
nosity. This difference indicates that the radio emission in some 
IsOQs might be mainly powered by the AGN. However, com- 
pared to HzRGs, IsOQs are less radio luminous, and they show 
higher MIR/radio luminosity ratios (see Table Q]). The presence 
of AGN-driven intermediate radio activity in some IsOQs has 
also been found by a study on a sample of MIR selected sources 
that include some of t he IsOQs in our sa mple and others at 
lower MIR luminosity (Saii na et alj|20()7h). This study claims 
that some IsOQs are at the beginning of a radio-loud phase. In 
conclusion, we find a variety of radio properties and no signif- 
icant correlation with the L(MIR)/L(NIR) ratio in the selected 
AGN samples. 



2.4. Star formation rates 

Recent evolutionary models also postul ate an interplay between 
AGN activity and SF (e.g. ISilkl 120051) . Interestingly, in spite 
of the difference in the stellar content, AGN power, and radio- 
activity in these systems, all SMG/AGNs, many HzRGs, and 
some IsOQs are characterized by high SFRs. The SFRs of 8 
IsOQs, obtained using MIP S FIR (70, and 160 u rn) data, range 
from 600 to 3,000 Moyr" 1 dPolletta et al.l l2008). Ty pical SFRs 
for S MG/AGNs range from 200 to 2,000 M^r" 1 (jPope et al.l 
2006). HzRGs are also characterized by high SFRs, up to a few 
1,000 M yr~' dReuland et al.ll2004t lArchibald et al.ll200ll) . The 
similarity in SFRs can be interpreted as lack of evidence for ra- 
dio activity or for AGN radiative power as regulators of SFR. 
This suggests either that AGN feedback does not take place as 
radio activity or AGN radiation, that its effects on the SFR are 
not observable yet, or that the feedback process does not act 
directly on the SFR. A similar result has recently been found 
by studying a sample of X-ray selected AGNs whose UV-NIR 
emiss ion is dominated by stellar light (lAlonso-Herrero et al.l 
2008). In this study, no evidence of either suppressed or elevated 
SF is found in galaxies with AGN activity compared to galax- 
ies of similar stellar masses and redshifts, suggesting that AGN 
activity does not affect the SFR of their host-galaxies. 
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Fig. 2. Median SEDs of IsOQs (red line in panel a), HzRGs (blue line in panel b), and SMG/AGNs (green line in panel c) normalized 
at their median MIR luminosity. The shaded areas are obtained from the 50% of all points closer to the median SED. The light-blue 
dotted line represents an extincted QSO template with Ay as indicated. The magenta dashed line represents the starburst (SB) galaxy 
M 82 with L(H) equal to the median stellar L(H) of each group, as annotated. The purple curves re present the sum of the extincted 
QSO template and of the starburst template. The Galactic center extinction curve (IChiar et al.ll2006h was applied to redden the QSO 
template. 



3. Discussion and conclusions 

We analyzed the stellar/ AGN luminosity ratio, the level of radio 
activity, and the SFR in three classes of AGNs at z~2, IsOQs, 
HzRGs, and SMG/AGNs, to investigate the existence of a link 
between AGN-driven radio activity, the buildup of stellar mass, 
and BH growth. We find that these AGNs are characterized by 
a broad range of stellar/AGN (L(NIR)/L(MIR)) luminosity ra- 
tios. On average, SMG/AGNs and HzRGs have more massive 
hosts than the most luminous, obscured QSOs at similar red- 
shifts. Such a broad range can be explained either by different 
Eddington ratios or by different Mb//-M» relations. 

Here, we have considered the possibility of an evolution- 
ary link among these various AGNs, but we do not rule out the 
possibility that the various types of AGNs are the result of dif- 
ferent evolutionary histories or physical states. Assuming that 
the stellar and BH mass can only grow, any link among these 
three classes of AGNs implies either that IsOQs will become 
less AGN-luminous and their hosts more massive, as observed 
in SMG/AGNs, or that SMG/AGNs will re semble IsOQs, as 
their AGN becomes m ore powerful (see e.g. iBorvs et aT]|2005t 
lAlexander et al .12005 al) . This last hypothesis is not supported by 
the relatively low stellar masses measured in IsOQs. However, 
it is possible that our IsOQ sample is biased towards systems 
with particularly low-mass hosts, as they are selected on the ba- 
sis of red optical-IR colors, we cannot rule out the latter sce- 
nario. On the other hand, the former scenario, evolution from 
IsOQs to SMG/AGNs, is quite plausible and supported by other 
observations. This scenario would imply that the bulk of stellar 
mass is still being assembled in IsOQs, while the bulge is at a 
more advanced stage in HzRGs and SMG/AGNs. HzRGs and 
SMG/AGNs might thus represent a later phase in the AGN evo- 
lution than are IsOQs. According to this scenario, fully grown 
BHs would already be in place in high-z galaxies before the 
bulk of stellar mass was assembled. This scenario is also con- 
sistent with observations of unobscured QSOs at very high-z 
(z >5) in which a fully grown BH and large reservoirs of molec- 
ular gas are present, but w here spheroidal stars have not yet 
formed dWalter et alj 120041) . Another supportive piece of evi- 
dence for this scenario is given by the younger stellar popu- 
lations found in unobscured AGNs compared to radio galax- 



Table 1. Radio properties of selected AGNs 



Source name 


z 


T a 
M.4Gi?9_ 


Log( f ) 






(W Hz" 1 ) 




IsOQ/LH_A4 


2.54 


25.04+0.02 


5.27 


IsOQ/LH_A5 


1.94 


24.54 


±0.02 


4.97 


IsOQ/LH_A6 


2.20 


24.42+0.03 


5.09 


IsOQ/LH_A8 


2.42 


26.71 


+0.02 


3.22 


IsOQ/N2_09 


1.98 


25.69 


+0.03 


4.15 


MIPS22204 


2.08 


25.68 


+0.02 


4.30 


Median(IsOQ£, 

Range(IsOQ)'' 


2.14+0.19 


25.68 


+0.68 


4.97+0.61 


1.94-2.54 


24.42- 


-26.71 


3.22-5.27 


Median(HzRG) 


1.90±0.53 


27.89 


+0.48 


1.11±0.37 


Range(HzRG) 


1.08-3.09 


26.55- 


28.53 


0.60-2.30 


Median(SMG/AGN) 


2.20±0.43 


24.31 


+0.28 


4.25±0.39 


Range(SMG/AGN) 


0.98-2.58 


23.25- 


24.63 


3.71-5.18 



a The radio lum inosities are from the deep VLA survey of the 
Lockman Hole dPolletta et al] 12006, Owen et al., in prep.) for 
LH_A4, LH.A5, LH_A6, and LH.A8, fro m the FIRST survey for 
N2.09, and from the VLA-FLS survey dSaiina et al.1 l2007bl) for 
MIPS22 204. The radio lumin ositi es for the HzRGs and SMG/AGNs 
are from Seymour et al. (2007) and Pope et al. ( 2006), respectively. 

* Includes only the 6 IsOQs with available radio data. The median 
redshift for the whole IsOQ sample is 2.23+0.31, and the redshift range 
is 1.29-2.96. 

tes dKauffmann et al.1 12008). However, this result has been ob- 
tained by studying low-z samples where we can expect different 
processes to dominate and evolutionary sequences to take place. 
Better defined samples and more accurate stellar and BH masses 
are needed to distinguish the two possible scenarios or rule out 
any evolutionary link among the three selected AGN samples. 

Finally, we want to underline the risks of assuming 
Eddington-limited AGN emission or the local Mj^-M, relation 
in the study of AGN and galaxy evolution. 
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